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Blending two discrete integrability criteria: 
singularity confinement and algebraic entropy 

S. Lafortune, A. Ramani, B. Grammaticos, 
Y. Ohta and K.M. Tamizhmani 



Abstract. We confront two integrability criteria for rational mappings. The 
first is the singularity confinement based on the requirement that every singu- 
larity, spontaneously appearing during the iteration of a mapping, disappear 
after some steps. The second recently proposed is the algebraic entropy cri- 
terion associated to the growth of the degree of the iterates. The algebraic 
entropy results confirm the previous findings of singularity confinement on 
discrete Painleve equations. The degree-growth methods are also applied to 
linearisable systems. The result is that systems integrable through linearisa- 
tion have a slower growth than systems integrable through isospectral methods. 
This may provide a valuable detector of not just integrability but also of the 
precise integration method. We propose an extension of the Gambier mapping 
in N dimensions. Finally a dual strategy for the investigation of the integra- 
bility of discrete systems is proposed based on both singularity confinement 
and the low growth requirement. 



1. Introduction 

The detection of discrete integrability has for the past few years relied on the 
criterion of singularity confinement [1]. This criterion is based on the observation 
that for integrable discrete systems any singularity that appears spontaneously 
disappears after some iterations. Let us give one example. In the mapping 

a 1 

x n+1 +x n - 1 = h — , (1.1) 

Xn X n 

where a is constant, we assume that, at some iteration, x n -\ vanishes while x n - 2 is 
finite. This has as consequence a diverging x n , and iterating once more a vanishing 
x n+ \. If we try to compute x n+2 we obtain an indeterminate form, oo — oo. In order 
to lift the indeterminacy we introduce a small parameter e and assume that x n -\ = 
e. Iterating we obtain successively, expanding in e, x n = 1/e 2 + a/e — x n - 2 + 0(e), 
x n +i = -e+ ae 2 + C(e 3 ), x n+2 = x n _ 2 + 2(ax n - 2 + l)e + C(e 2 ). By taking 
the appropriate limit e — > we find that x n + 2 is finite, and moreover contains 
the information on the initial data i.e. x n _ 2 . The subsequent indeed 
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finite and the singularity is confined to the sequence {0,oo 2 ,0}, the square in the 
infinite sign being a reminder of the 1/e 2 . Singularity confinement is tailored for 
rational mappings although this should not be considered as an absolute restriction. 
What singularity confinement cannot treat are mappings which do not possess 
singularities. For example a mapping of the form x n+ i = Xx n (l — x n ) lies beyond 
the reach of the confinement method [2]. Our tacit conjecture has always been that 
polynomial mappings, like the logistic one, are not integrable, with only exception 
the linear one [3]. 

While confinement has turned out to be a most useful integrability detector, its 
application was based on the unwarranted assumption that its necessary character 
will be constraining enough to make it sufficient. Still, from the outset it was 
clear that confinement was not a sufficient criterion. Let us recall the example 
presented already in [4]. If we examine the singularity structure of the mapping 
x n+ i = f(x n ,n) where the / is rational in x n and analytic in n we find that all 
f's of the form / = J2k ( x +/3 k )"k > v k € Z, lead to confinement. Clearly 

not all these mappings can be integrable and, indeed, only the homographic one, 
x n +\ = oe+ - q-g , is. The non sufficient character of singularity confinement was put 
to a wider perspective by Hietarinta and Viallet [5] who verified that the mapping: 

1 / X 

x n+1 + X n -! = x n + — (1.2) 

X n 

is confining, the singularity sequence being {0, oo 2 , oo 2 , 0}, and at the same time 
exhibits chaotic behaviour. Moreover this is not an isolated example and one can 
easily construct mappings which are confining without being integrable. 

Clearly a more stringent integrability criterion was needed and the authors of 
[5] (see also [6]), have proposed one based on the ideas of Arnold and Veselov. 
According to Arnold [7] the complexity (in the case of mappings of the plane) is 
the number of intersection points of a fixed curve with the image of a second curve 
obtained under the mapping at hand. While the complexity grows exponentially 
with the iteration for generic mappings, it can be shown [8] to grow only polynomi- 
ally for a large class of integrable mappings. As Veselov points out, "integrability 
has an essential correlation with the weak growth of certain characteristics" . Thus 
the authors of [5] proposed to directly test the degree of the successive iterates and 
introduced the notion of algebraic entropy. It is defined as E — lim^oo log(d„)/n 
where d n is the degree of the nth iterate of some initial data under the action of 
the mapping. Since a generic nonintegrable mapping exhibits exponential degree- 
growth, a nonzero algebraic entropy indicates nonintegrability. Integrable mappings 
must have zero algebraic entropy, associated to slower-than-exponential (typically 
polynomial) degree-growth. It is expected that the requirement of zero algebraic 
entropy is strong enough to be a sufficient integrability condition. 

In order to make the ideas clearer we present here our method for the study 
of the degree-growth. It presents a few differences with respect to the one of [5] 
which, we feel, make its practical implementation simpler. We start by introducing 
homogeneous variables through the appropriate choice of initial data. Typically, in 
a three-point mapping we choose to introduce xq = r, x\ = p/q and thus the weight 
of r is while that of p, q are taken equal to 1. We then compute the homogeneous 
degree in p, q in the denominator and numerator of x n at every iteration. Let us 
follow the first few iterations of mapping (1.1) in order to understand the mechanism 



BLENDING TWO DISCRETE INTEGR ABILITY CRITERIA 



3 



of the degree-growth. We find 



q + apq — rp pP^ 



X2 = 5 ) x 3 



p 2 ' q(q 2 + apq — rp 2 ) 2 ' 

_ (g 2 + apg-rp 2 )P 6 _ P 4 P 9 

X4 ~ p2 ' ^5 — p2 

where the P^'s are homogeneous polynomials in p, q of degree k. (Remember that 
r is of zero homogeneous degree in our convention). A pattern becomes apparent. 
Whenever a new polynomial appears in the numerator of x n , its square will appear 
in the denominator of x n +i and it will appear one last time as a factor of the 
numerator of x n +2, after which it disappears due to factorisations. The singularities 
we are working with in the singularity confinement approach correspond to the zeros 
of any of these polynomials, which explains the pattern {0, oo 2 , 0}. The singularity 
confinement is intimately related to this factorisation which plays a crucial role in 
the algebraic entropy approach. If we calculate the degree of the iterates, we obtain: 
0, 1, 2, 5, 8, 13, 18, 25, 32, 41, ... Clearly the degree-growth is polynomial: we 
have dim — 2™ 2 and d 2m +i = 2m 2 + 2m + 1. (A remark is necessary at this point. 
In order to obtain a closed-form expression for the degrees of the iterates, we start 
by computing a sufficient number of them. Once the expression of the degree has 
been heuristically established we compute the next few ones and check that they 
agree with the analytical expression predicted) . Thus the algebraic entropy of this 
mapping is zero in agreement with its integrable character [9] . 

On the other hand, if we iterate mapping (1.2) we obtain the sequence 

p 3 — p 2 qr + q 3 p - 



X2 = n , £3 



p2g p 2 (p 3 — p 2 qr + q 3 ) 2 ' 

_ P-P22 _ {p 3 -p 2 qr + g 3 )P 5 8 

Xi ~ q(p 3 - p 2 qr + q 3 ) 2 P 8 2 ' ^ ~ qP s 2 P 22 2 

Again some factorisations and simplifications do occur which explain why this map- 
ping has the confinement property. However the degree of the new terms appearing 
at every iteration grows too rapidly and thus the simplifications cannot curb the 
exponential growth. The degree sequence now is 0, 1, 3, 8, 23, 61, 160, 421, . . . , i.e. 
the degrees obey the relation d n +i — 3d n + d n -\ = |(1 + j n+1 + j 2< - n+1 )) where j is 
a complex cubic root of 1 leading to an exponential growth with asymptotic ratio 
(3 + \/5)/2 and algebraic entropy log((3 + Vb)/2). The situation is even worse in 
the case of nonconfining mappings where no simplifications occur and the growth 
rate is maximal. 

Since a considerable volume of results on integrable discrete systems were ob- 
tained with the method of singularity confinement it is natural to question their 
integrability in the light of the findings concerning the non sufficiency of the sin- 
gularity confinement criterion. One of the aims of this review is to provide the 
confirmation of these results based on the more stringent criterion of algebraic 
entropy. In particular we shall examine several discrete Painleve equations that 
we have obtained over the years and obtain their growth properties. The main 
result of this analysis is that, in every case studied, the singularity confinement 
criterion when used for deautonomisation of an integrable mapping, turns out to 
be sufficient [10]. This means that an equation, the autonomous form of which is 
integrable and which has been deautonomised following the constraints provided 
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by singularity confinement, leads to exactly the same degree-growth as the initial, 
autonomous, one. Thus we expect the deautonomised equation to be as integrable 
as its autonomous limit. This finding is of capital importance since it validates the 
results previously obtained with the singularity confinement approach. Moreover 
it sets the frame for a new, composite, approach. Whenever we wish to obtain 
the integrable cases of an equation containing a certain number of parameters we 
can perform a first, exploratory, study using the confinement method. Once the 
confinement constraints have been used in order to limit the freedom, the algebraic 
entropy, or low-growth criterion, can be implemented in order to pin down the truly 
integrable cases. Examples of this dual approach can be found in [11]. 



2. Discrete Painleve equations 

Let us first recall what has always been our approach to the derivation of dis- 
crete Painleve equations. We start from an autonomous system the intcgrability of 
which has been independently established. In the case of discrete Painleve equa- 
tions, this system is the QRT mapping [9]: 

f (1) (xn) - (x n+ i + x n - 1 )f {2 \x n ) + x n+1 x n ^f^(x n ) = (2.1) 

When the /W's are quartic functions, satisfying specific constraints, the mapping 
(2.1) is integrable in terms of elliptic functions. Since the elliptic functions are the 
autonomous limits of the Painleve transcendents, the mapping (2.1) is the appro- 
priate starting point for the construction of the nonautonomous discrete systems 
which arc the analogues of the Painleve equations. The procedure we used, of- 
ten referred to as 'deautonomisation', consists in finding the dependence of the 
coefficients of the quartic polynomials appearing in (2.1) with respect to the inde- 
pendent variable n, which is compatible with the singularity confinement property 
Namely, the n-dependence is obtained by asking that the singularities are indeed 
confined. The reason why this procedure can be justified is the following. Since 
the autonomous starting point is integrable, it is expected that the growth of the 
degree of the iterates is polynomial. Now it turns out that the application of the 
singularity confinement deautonomisation corresponds to the requirement that the 
nonautonomous mappings lead to the same factorizations and subsequent simpli- 
fications and have precisely the same growth properties as the autonomous ones. 
These considerations will be made more transparent thanks to the examples we 
present in what follows. 

Let us start with a simple case. We consider the mapping: 



a 1 .„ 

x n+1 +x n -! = h— (2.2) 

x n x n 



where a depends on n. The singularity confinement result is that a must satisfy the 
conditions a n+ \ — 2a n + a„_i = i.e. a is linear in n. Assuming now that a is an 
arbitrary function of n we compute the iterates of (2.2). We obtain the sequence: 

q 2 + aipq - p 2 r pQ 4 

Xl = 7i , X 3 



p 2 q(q 2 + a\pq — p 2 r) 2 ' 

(q 2 + aipq - p 2 r)Q 7 pQaQ 



X\ — ^-.o j x§ — 



12 



pQl ' q(q 2 + mpq - p 2 r)Qj 
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where the Qfc's are homogeneous polynomials in q, r of degree k. The simplifications 
that do occur are insufficient to curb the asymptotic growth. As a matter of fact, 
if we follow a particular factor we can check that it keeps appearing either in the 
numerator or the denominator (where its degree is alternatively 1 and 2). This 
corresponds to the unconfincd singularity pattern {0,oo 2 ,0,oo,0,oo 2 ,0,oo,...}. 
The confinement condition a n+ i — 2a n + a„_i = is the condition for q to divide 
exactly Qi, for both q and r 2 + a\qr — pq 2 to divide exactly Q12, etc. Let us 
now turn to the computation of the degrees of x n . We obtain successively 0, 1, 
2, 5, 10, 21, 42, 85, .... The growth is exponential, the degrees behaving like 
rf2m-i = (2 2m — l)/3 and d,2m = 2d 2 m-i, a clear indication that the mapping 
is not integrable in general. Already at the fourth iteration the degrees differ in 
the autonomous and nonautonomous cases. Our approach consists in requiring 
that the degree in the nonautonomous case be identical to the one obtained in the 
autonomous one. If we implement the requirement that d± be 8 instead of 10 we 
find the condition a n+ i — 2a„ + a„_i = 0, i.e. precisely the one obtained through 
singularity confinement. Moreover, once this condition is satisfied, the subsequent 
degrees of the nonautonomous case coincide with that of the autonomous one. 
Thus this mapping, leading to polynomial growth, should be integrable, and, in 
fact, it is. As we have shown in [9], where we presented its Lax pair, equation 
(2.2) with a n = an + {3 is a discrete form of the Painlcvc I equation. In the 
examples that follow, we shall show that in all cases the nonautonomous form of 
an integrable mapping obtained through singularity confinement leads to exactly 
the same degrees of the iterates as the autonomous one. 

We now turn to what is known as the "standard" discrete Painleve equations 
[12] and compare the results of singularity confinement to those of the algebraic 
entropy approach. We start with d-Pi in the form: 



The degrees of the iterates of the autonomous mapping are 0, 1, 2, 3, 6, 9, 12, 
17, 22, . . . , i.e. a quadratic growth with <i3 m +fc=3m 2 + (2m + l)fc, for k =0,1,2, 
while those of the generic nonautonomous one are 0, 1, 2, 3, 6, 11, ... Requiring 
two extra factorisations at that level (so as to bring d 5 down to 9) we find the 
following conditions a n+ \ = a n , so a n must be a constant, and & n +2 — b n +i — 
b n + b n -i=0, i.e. b n is of the form b n = an + j3 + -f(— 1)" which are exactly the 
result of singularity confinement. Implementing these conditions we find that the 
autonomous and nonautonomous mappings have the same (polynomial) growth 
[10]. Both are integrable, the Lax pair of the nonautonomous one, namely d-Pi 
having been given in [13,14,15]. 

For the discrete Pu equation we have 



The degrees of the iterates in the autonomous case are d n =0, 1, 2, 4, 6, 9, 12, 16, 20, 
. . . , (i.e. d2m-i~fn 2 , d 2 m=n^ 2 + m) while in the generic nonautonomous case we 
find the first discrepancy for d^ which is now 8. To bring it down to 6 we find two 
conditions, a„ +1 — 2a„ + a„_!=0 and b n+1 = 6„_i. This means that a is linear in 
n and b is an even/odd constant, as predicted by singularity confinement. Once we 
implement these constraints, the degrees of the nonautonomous and autonomous 




(2.3) 



ax n + b 



(2.4) 
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cases coincide. The Lax pair of equation (2.4) in the nonautonomous form, i.e. 
d-Pu, has been presented in [13,16,17]. 

The q-Pui equation was obtained from the deautonomisation of the mapping: 

_ (x n - a){x n - b) 
Xn+lXn - 1 ~ (l-cx n )(l-x n /c) (2 ' 5) 

In the autonomous case we obtain the degrees (i n =0, 1, 2, 5, 8, 13, 18, . . . , just 
like for equation (2.2), while in the generic nonautonomous case we have 0, 1, 2, 
5, 12, ... For d± to be 8 instead of 12, one needs four factors to cancel out. The 
conditions are c„ + i = c„_i and a„ + i&„_i = a„_i6„ + i = a n b n . Thus c is a constant 
up to an even/odd dependence, while a and b are proportional to A" for some A, 
with an extra even/odd dependence, just as predicted by singularity confinement 
in [12]. The Lax pair for q-Pm has been presented in [13,18]. 

For the remaining three discrete Painleve equations the Lax pairs are not known 
yet. It is thus important to have one more check of their integrability provided by 
the algebraic entropy approach. We start with d-Piy in the form: 

(x n+1 +x n )(x n . 1+ x n ) = {X }~f} { f?~P (2.6) 

yx n -|- z n ) c 

where a, b and c are constants. If z n is constant we obtain for the degrees of the 
successive iterates d n —0, 1, 3, 6, 11, 17, 24, . . . The general expression of the growth 
is d n =6m 2 if n = 3m, rf„=6m 2 +4m+l if n = 3m +1 and d n =6m 2 + 8m + 3 if n = 
3m + 2. This polynomial (quadratic) growth is expected since in the autonomous 
case this equation is integrable, its solution being given in terms of elliptic functions. 
For a generic z n we obtain the sequence d n —0, 1, 3, 6, 13, . . . The condition for 
the extra factorizations to occur in the last case, bringing down the degree d^ to 
11, is for z to be linear in n. We can check that the subsequent degrees coincide 
with those of the autonomous case. 

For the discrete Painleve V equation we start from: 

(x n+1 x n - l)(x n . lXn - 1) = (zl + ^ + l)(xl + bx + 1) (2 ?) 

yL ZnCX n )yl ZndX n J 

where a, 6, c and d are constants. If moreover z is also a constant, we obtain exactly 
the same sequence of degrees d n =0, 1, 3, 6, 11, 17, 24, . . . , as in the d-Piy case. 
Again, this polynomial (quadratic) growth is expected since this mapping is also 
integrable in terms of elliptic functions. For the generic nonautonomous case we 
again find the sequence d„=0, 1, 3, 6, 13, ... Once more we require a factorization 
bringing down cLt to 11. It turns out that this entails a z which is exponential in 
n, which then generates the same sequence of degrees as the autonomous case. In 
both the d-Piy and q-Pv cases we find the n-dependence already obtained through 
singularity confinement. Since this results to a vanishing algebraic entropy we 
expect both equations to be integrable. 

The final system we shall study is the one related to the discrete Pyi equation: 

(x n -\-ix n z n ^\z r ^){x n —\x n z n _iz n ) {x n ~\~ dz n x n + z n ){x ri -\- bz n x n -\- z n ) 

1) (xl + cx n + l)(xl + dx n + 1) 

(2.8) 

where a, b, c and d are constants. In fact the generic symmetric QRT mapping 
can be brought to the autonomous (z n constant) form of equation (2.8) through 
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the appropriate home-graphic transformation. In the autonomous case, we obtain 
the degree sequence d n =0, 1, 4, 9, 16, 25, ... , i.e. d n =n 2 . Since mapping (2.8) 
is rather complicated we cannot investigate its full freedom. Still we were able to 
perform two interesting calculations. First, assume that in the rhs instead of the 
function z n a different function ( n appears. In this case the degrees grow like 0,1, 
5, . . . , and the condition to have d 2 =4: instead of 5 is z n+ iz n -iz 2 = Assuming 
this is true, we compute the degree c?3 of the next iterate and find ds=13 instead of 
9. To bring down c?3 to the value 9 we need z 2 = ( 2 , which up to a redefinition of 
a and b means z n = ( n . This implies z„ + iz„_i = z 2 , and z n is thus an exponential 
function of n, z n =z \ n (which is in agreement with the results of [19]). Then 
a quartic factor drops out and c? 3 is just 9. One can then check that the next 
degree is 16, just as in the autonomous case. Thus the q-Pvi equation leads to the 
same growth as the generic symmetric QRT mapping and is thus expected to be 
integrable. As a matter of fact we were able to show that the generic asymmetric 
QRT mapping leads to the same growth d n —n 2 as the symmetric one. This is not 
surprising, given the integrability of this mapping. What is interesting is that the 
growth of the generic symmetric and asymmetric QRT mappings are the same. 
Thus d n =n 2 is the maximal growth one can obtain for the QRT mapping in the 
homogeneous variables we are using. As a matter of fact we have also checked that 
the asymmetric nonautonomous q-Pyi equation, introduced in [19] led to exactly 
the same degree-growth d n =n 2 . 

The singularity confinement results have been confronted to the algebraic en- 
tropy approach for several other discrete Painleve equations. In every case examined 
the deautonomisation obtained has turned out to be the right one, it was the con- 
dition for the degree-growth to be identical to the one of the autonomous system. 
Thus despite the non sufficiency of the singularity confinement the integrability 
predictions for discrete Painleve equations based on this criterion are confirmed. 

3. Linearisable equations 

In this section we shall examine this particular class of mappings which are 
linearisable and study their growth properties. Most of these systems were ob- 
tained using the singularity confinement criterion and thus a study of the growth 
of the degree of the iterates would be an interesting complementary information. 
Moreover, as we will show, the linearisable systems do possess particular growth 
properties which set them apart from the other integrable discrete systems. 

The first mapping we are going to treat is a two-point mapping of the form 
x n+ \ = f(x n ,n) where / is rational in x n and analytical in n. As explained in 
the introduction for all f's of the form J2k ( x +p k )"k tnc singularity confinement 
requirement is satisfied. However only the discrete Riccati, x n +i = a + - x +/3 , is 
expected to be integrable. Our argument in [20], for the rejection of these confining 
but nonintegrable cases, was based on the proliferation of the preimages of a given 
point. If we solve the mapping for x n in terms of x n+ i we do not find a uniquely 
defined x n and, iterating, the number of x n _k grows exponentially. In what fol- 
lows we shall analyse this two-point mapping in the light of the algebraic entropy 
approach. We start from the simplest case which we expect to be nonintegrable, 



(3.1) 
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The initial condition we are going to iterate is xo = p/q and the degree we calculate 
is the homogeneous degree in p and q of the numerator (or the denominator) of the 
iterate. We obtain readily the following degree sequence d n — 1, 2, 4, 8, 16, . . . i.e. 
d n = 2™. Thus the algebraic entropy of the mapping is log(2) > 0, an indication 
that the mapping cannot be integrable. Now we ask how can one curb the growth 
and make it nonexponential. It turns out that the only possibilities are = 
or (3 = 7. In cither case mapping (3.1) becomes a homography. The degree in 
this case is simply d n = 1 for all n. This is an interesting result, clearly due to 
the fact that the homographic mapping is linearisable through a simple Cole-Hopf 
transformation. 

The second mapping we shall examine is one due to Bellon and collaborators 

[21] 



X n +l 



Vn+1 = 



Vn {%n Un ) 



(3.2) 



The degree-growth in this case is studied starting from x = r, y = p/q and again 
we calculate the homogeneous degree of the iterate in p and q, i.e. we set the degree 
of r to zero. (Other choices could have been possible but the conclusion would not 
depend on these details.) We obtain the degrees d Xn = 0,2,2,4,4,6,6,... and 
d Vn = 1,1,3,3,5,5,... i.e. a linear degree-growth. This is in perfect agreement 
with the integrable character of the mapping. As was shown in [22] it does satisfy 
the unique preimage requirement and possesses a constant of motion k = 
the use of which reduces it to a homographic mapping for x n or y n . 

The third mapping we are going to study is the one proposed in [20] 



XnVn 



X n +1 



Xn {,X n Vn 



Xn Vn (3.3) 

(X n - y n ){x n ~ Vn ~ a) 

Vn+1 — Y~Z 

x n Vn 

where a was taken constant. We start by assuming that a is an arbitrary function 
of n and compute the growth of the degree. We find d Xn = 0, 1, 2, 3, 4, 5, 6, 7, 8, . . . 
and d Vn = 1, 2, 3, 4, 5, 6, 7, 8, 9, . . . i.e. again a linear growth. This is an indication 
that (3.3) is integrable for arbitrary a n and indeed it is. Dividing the two equations 
we obtain y n+1 /x n+ i = 1 - y n /x n i.e. y n /x n = 1/2 + fc(-l)" whereupon (3.3) is 
reduced to a homographic mapping for x. Thus in this case the degree-growth has 
successfully predicted integrability. 

A picture starts emerging at this point. While in our study of discrete Painleve 
equations and the QRT mapping we found quadratic growth of the degree of the 
iterate, linearisable second-order mappings seem to lead to slower growth. In order 
to investigate this property in detail we shall analyse the three-point mapping 
we have studied in [4,23] from the point of view of integrability in general and 
linearisability in particular. The generic mapping studied in [23] was one trilinear 
in x n , x n+ i , x n -\. Several cases were considered. Our starting point is the mapping, 

x n+1 x n x n - 1 +(3x n x n+1 +(r]x n+1 x n - 1 +ix n x n - 1 +0ix n + rix n - 1 +(x n+1 + 1 = 0. 

(3.4) 
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We start with the initial conditions x = r, x\ = p/q and compute the homoge- 
neous degree in p, q at every n. We find d n = 0, 1, 1, 2, 3, 5, 8, 13, . . . i.e. a Fibonacci 
sequence d n+ i = d n + d n -\ leading to exponential growth of d n with asymptotic 
ratio 1+ 2 V ^ ■ Thus mapping (3.4) is not expected to be integrable in general. How- 
ever, as shown in [23] integrable subcases do exist. We start by requiring that the 
degree-growth be less rapid and as a drastic decrease in the degree we demand that 
rf 3 = 1 instead of 2. We find that this is possible when either /3 = £ = in which 
case the mapping reduces to: 

n 1 

x n +i = -7 (3-5) 

Xn X n X n —± 

or 7 = i] = 0, giving a mapping identical to (3.5) after x — > 1/x. In this case the 
degree is d n = 1 for n > 0. Equation (3.5) is the generic projective three-point 
mapping, written in canonical form. Its linearisation can be obtained [23] in terms 
of a system of three linear equations, a fact which explains the constancy of the 
degree. 

Non generic subcases of (3.4) some of which are integrable do exist. They have 
been studied in [23,24]. 

Linearisability through the reduction to a projective system is not the only 
possibility. Other possibilities do exist. Let us start by considering the generic 
three-point mapping that can be considered as the discrete derivative of a (discrete) 
Riccati equation. Let us start from the general homographic mapping which we 
can write as 

Ax n x n+ i + Bx n + Cx n+ i + D = (3.6) 
where A, B, C, D are linear in some constant quantity k. In order to take the 
discrete derivative we extract the constant k and rewrite (3.6) as: 

QLXnXn+l 

k= -. (3.7) 

eXnXn+l + Qx n + T}X n Up + 6 

Using the fact that n is a constant, it is now easy to obtain the discrete derivative by 
downshifting (3.7) and subtracting it form (3.7) above. Instead of examining this 
most general case we concentrate on the forms proposed in [25]. They correspond 
to the reduction of (3.7) to the two cases: 

b , , 

k = x n+ i + a-\ (3.8) 

Xn 

K= x n+1 (x n + a) 

X n +b 

Next we compute the discrete derivatives of (3.8) and (3.9). We find: 

bn b n — i 

x n +i = x n + a n -i - a n 1 (3.10) 

and 

%n— 1 H~ a n— 1 %n H~ /o 1 1 \ 

X n +l=X n — . (3.11) 

Xn i Q>n X n —1 T On—I 

The study of the degree of growth of (3.10) and (3.11) can be performed in a 
straightforward way. For both mappings we find the sequence d n = 0, 1, 2, 3, 4, 5,6,.. 
i.e. a linear growth just as in the cases of mappings (3.2) and (3.3). If we substitute 
b n -i by c n _i in the last term of the rhs of (3.10) or the denominator of (3.11) we 
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find d n = 0, 1, 2, 4, 8, 16, . . . i.e. d n — 2 n for n > unless c = b. Investigating all 
the possible ways to curb the growth we find for both (3.10) and (3.11) that c = 
is also a possibility to bring (i 3 down to 3. However a detailed analysis of this case 
shows that for c = we have d n = 0, 1, 2, 3, 5, 8, 13, 21, . . . i.e. a Fibonacci sequence 
with slower, but still exponential, growth (i.e. ratio 1+ 2 V ^ instead of 2). 

A bonus study: the N- dimensional Gambier mapping 

While most of the results included in this paper were presented in previous 
publications of ours, in this section we shall present a study which appear for the 
first time in these proceedings. It concerns the extension of the Gambier mapping 
to N dimensions. 

The two-dimensional Gambier mapping was introduced in [26,27] as a discreti- 
sation of the second order differential equation discovered by Gambier, in his study 
of second order ODE's having the Painleve property. The Gambier equation is in 
fact a system of two Riccati's in cascade. The latter means that the system consists 
in one Riccati for one variable and a second one (for the second variable) with co- 
efficients depending (linearly) on the variable of the first one. The discretisation of 
the Gambier system is straightforward: just substitute homographic mapping, i.e. 
discrete Riccati, in the previous sentence. Thus the general second-order Gambier 
mapping is given by: 

ay n + b 

= ^Td' (3 ' 12a) 

Xn+1 = f Vn l% Xn l^ n Z h \ (3.12&) 
\3Vn + k)x n + (ly n + to) 

where a, b, . . . , to are functions of n. In [27] we have studied this mapping in 

detail from the point of view of the singularity structure. This was motivated by 

the fact that we aimed at being able to express the solution as an infinite product 

of matrices, even across singularities. On the other hand if we are not interested 

in this fine point, the linearisability of (3.12) can be obtained through a Cole-Hopf 

transformation for each variable. 

The study of the degree-growth of (3.12) is straightforward. We start from 
xq = r, yo = p/q and compute the homogeneous in p, q degree of (4.1a) and (4.1b). 
Since (3.12a) is a Riccati its degree does not grow i.e. we have d Vn = 1. Given 
the structure of (3.12b) we have d Xn+1 — d Xn + d yn and thus d Xn — n. What 
is interesting here is that the Gambier mapping exhibits a linear degree-growth 
independently of the precise values of a, b, . . . , to. The fact that it can be reduced 
to Riccati's in cascade is enough to guarantee its integrability. 

The generalisation of (3.12) to N + 1 dimensions is straightforward. We find 
the system: 

(o) _ a^x^ + M°> 
Xn+1 cWsW+dW 

x (i) _ (e^xt 1} + /Wjgjg + (g^xt 1} + ftW) 

Xn+1 (j^xt 1} + k^)x { : ] + {mx { t 1] + mW) ' 
The study of the degree-growth of system (3.13) can be performed along the 
N = 1 case. From (3.13) we have the recurrence: G?a4+i = dx^ + dx^n ^ • We 
obtain formally dx^ x = J2p=o dxp 1-1 ^ . Thus starting from dx„ ^ = 1 we find 



(3.13) 

l,...,N. 
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dx n cx n N . Thus the (N + l)-dimensional Gambier mapping has polynomial 
growth for every dimension. The linearisation is obtained through a Cole-Hopf 
transformation. 

Given the structure of (3.13) it is clear that we can solve successively for each 
variable and express finally (3.13) as a single (N + 2)-point mapping. This leads us 
to another integrable discretisation of linearisable mappings with cascade structure. 
Let us consider the second order system 



where the f^s are polynomial in y n . From the structure of (3.14) it is clear that it 
can be linearised, independently of the precise form of the fi . Similarly the degree- 
growth of x n is always linear. On the other hand it is not, in general, possible to 
express (3.14) as a single three-point mapping. Extension of (3.14) to TV-dimensions 
can be obtained along the lines of the iV-dimensional discrete Gambier system. 

In this section we have studied the growth properties of various linearisable 
systems identified through the singularity confinement criterion. In every case the 
singularity confinement results were confirmed. Moreover it turned out that the 
linearisable systems lead to slower growth than systems which are integrable by 
other methods. This property could be used for the classification of integrable 
systems and be a valuable indication as to the precise method of their integration. 



In this work we have presented a comparative review of results obtained with 
the singularity confinement and the algebraic entropy methods. While the for- 
mer approach is not a sufficient criterion of integrability, the confirmation of its 
results (through the second, more stringent, criterion) in the domain of discrete 
Painleve equations leads to interesting new insights. In every case examined, the 
singularity confinement, when used for the deautonomisation of an integrable au- 
tonomous mapping turned out to give sufficient constraints for the degree-growth 
to be non-exponential. Thus we can propose a new strategy for the detection of 
integrable discrete systems. Given a mapping which contains several parameters 
the singularity confinement necessary criterion can be used in order to screen it for 
possible integrable discrete cases. Once the research domain is reduced the growth 
properties can be studied leading to better integrability candidates. 

Another interesting result of our studies concerns the degree-growth of lin- 
earisable systems. We found that, while for second-order mappings the generic 
integrable case is associated to quadratic growth, the linearisable mappings lead to 
zero or linear growth. The growth exponent is of course a property which depends 
on the dimension. We surmise that the generic integrable iVth order mapping will 
lead to growth n N , while the linearisable mappings of the same dimension will lead 
to slower growth. Our study of the iVth order Gambier system shows that the 
growth is n"^ 1 , and we expect the the projective system of order N to lead to zero 
growth. Thus the detailed study of the growth properties can become a precious 
indication as to the precise method of integration of a given discrete system. 



X n +1 



Vn+1 



ay n + b 
cy n + d' 

fi(y n )x n - Myn) 
h{y n ) - h{y 



(3.14) 



4. Conclusion 
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